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ABSTRACT 



We present mid-IR interferometric observations of six type 1 AGNs at multiple baseline lengths ranging from 27 m to 130 m, 
reaching high angular resolutions up to A/B ~ 0.02 arcseconds. For two of the targets, we have simultaneous near-IR interferometric 
measurements as well, taken within a week. We find that all the objects are partially resolved at long baselines in these IR wavelengths. 
The multiple-baseline data directly probe the radial distribution of the material on sub-pc scales. We show that for our sample, which 
is small but spans over ~2.5 orders of magnitudes in the UV/optical luminosity L of the central engine, the radial distribution clearly 
and systematically changes with luminosity. 

The brightness distribution at a given mid-IR wavelength seems to be rather well described by a power law, which makes a simple 
Gaussian or ring size estimation quite inadequate. In this case, a half-light radius R112 can be used as a representative size. We show 
that the higher luminosity objects become more compact in normalized half-light radii Ri/2/Rm in the mid-IR, where is the dust 
sublimation radius empirically given by the L''~ fit of the near-IR reverberation radii. This means that, contrary to previous studies, 
the physical mid-IR emission size (e.g. in pc) is not proportional to L''^, but increases with L much more slowly. With our current 
datasets, we find that Rt/2 iP-^^^o.os g g ^j^^^ ^j^^j ^^^^ nearly constant at 13 fim. 

The derived size information also seems to correlate with the properties of the total flux spectrum, in particular the smaller /?i/2/^in 
objects having bluer mid-IR spectral shape. We use a power-law temperature/density gradient model as a reference, and infer that 
the radial surface density distribution of the heated dust grains at a radius r changes from a steep structure in high luminosity 
objects to a shallower ~/' structure in those of lower luminosity. The inward dust temperature distribution does not seem to smoothly 
reach the sublimation temperature - on the innermost scale of ~Ri„, a relatively low temperature core seems to co-exist with a slightly 
distinct brightness concentration emitting roughly at the sublimation temperature. 

Key words. Galaxies: active. Galaxies: Seyfert, Infrared: galaxies. Techniques: interferometric 



1. Introduction 

A key step in understanding the accretion process in active galac- 
tic nuclei (AGN) is to determine the role and structure of the pu- 
tative obscuring torus. The obscuration here essentially refers to 
the radiative absorption by dust grains in the torus, which hides 
the UV/optical emission of the central engine and the broad line 
emission (see Antonuccill 19931 for a review). In type 2 objects, 
which do not display these spectral components, our line of sight 
towards the nucleus is thought to be relatively close to the sys- 
tem's equatorial direction along which the torus obscuration be- 
comes effective. Type 1 objects are believed to be rather face-on 
and permit us to directly observe the central components un- 
obscured. The torus dust grains thermally re-emit the absorbed 
energy, and the innermost dust sublimates at ~1500 K, setting 
the inner boundary of the dust distribution. This dust sublima- 
tion region thus radiates mostly in the near-infrared (near-IR; 
wavelength A ~ a few yum), with outer dust grains emitting at 
lower temperatures in the mid-IR (~10 fim). Without much spa- 
tial information, the torus structure has long been inferred from 
the observed spectral energy distributions (SEDs) mainly in the 
near- and mid-IR, as well as theoretical arguments. 

The SEDs were first modeled by the tori with smooth 
densi ty distributions (Pier & Rrolik 1993; Granato & Danese 
Il994t lEf stathiou & Rowan-RobinsonI 119951) . However, the 
case for the material being in a clumpy medium has theoret- 



ically been argued from an early stage jKrolik & Begelmaiij 
Il988l) . Accordingly, there have been many eff orts to develop 
clumpy torus models in the past decade (iNenkova et aL 
2002; Dullemond & van Bemmel 2005; Honig et al. 
2006; Schartmann et al. i2008c iNenkova et all i2008ar 
Honig & Kishimoto 2010), and the torus structure has 
been inferred in the contex t of these clumpy toru s mod- 
Ramo s Almeida et al.' '2009t_ 
2010; Ramos Almeida et al 



2008bt 



els (INenkova et a 

M oretal.1 120091; iHonigetal _ 

i201ll Alonso-Herrero et al. n2oTir However, we are finally 
beginning to have far more direct information with the long- 
baseline interferometry in both the near- and mid-IR by spatially 
resolving the inner structure. 



1.1. Early Type 2 investigations 

Detailed interferometric studies in the mid-IR have targeted a 
few bright typ e 2 AG Ns, with the fi rst ones being NGC1068 
and Circinus (IJaffe et al.l[2004t iTristram et al J 120071) . using the 
mid-infrared interferometric instrument (MIDI) at the Very 
Large Telescope Interferometer ( VLTI). These have achieved the 
first determination of the overall size and projected orientation 
of the mid-IR emission region with two-temperature elliptical 
Gaussian models. The AGN NGC1068 was also successfully 
observed with the near-IR (2.2 yum) long-baseline interferomet- 
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ric instrument VINCI at VLTI (IWittkowski et alJl2004 . and fol- 
lowed up in the mid-IR with more detailed MIDI observations 
(Raban et al. 2009). One hard aspect is that, since the nucleus 
is heavily obscured in type 2 objects, one needs to disentan- 
gle the effects of radiative transfer in order to study the intrinsic 
structure and distribution of the material. Another type 2 object, 
Cen A, was the first radio-loud object to be studied with the IR 
long -baseline interferometry ( Meisenheimer et al. 2007), which 
was followed by more detailed observations bv iBurtscher et alj 
(I2OIO). Interestingly, the corresponding mid-IR source is found 
to be dominated by an unresolved source, interpreted to be a syn- 
chrotron core, which might, however, also complicate the study 
of the structure. 



1.2. Type 1 study in the near-IR 

Studies of type 1 objects are largely unaffected by the incli- 
nation/obscuration. We can directly scrutinize the inner intrin- 
sic structure in bo th the near-IR and mid-IR. In the near-IR, 
ISwain et al ] (I2003h made the first measurement of the K-band 
(2.2 jjm) visibility of the brightest type 1 AGN NGC4151 with 
the Keck interferometer (KI), where they obtained quite high 
visibility (squared visibili ty ~ 0.8). T his resu lt was re- 
cently confirmed by Kishi moto et al. (l2009al) a nd al so fPott et al.l 
(l2010l) . Furthermore, Kishi moto etal] ( 2009a. 2 01 ih found simi- 
larly high visibilities for in total eight type 1 AGNs. They argued 
that this is an indication of the marginal and partial resolution 
of the dust sublimation region in these type 1 AGNs. Firstly, 
they marginally detected a visibility decrease over the baseline 
length for NGC4151, which would be robust evidence for a par- 
tially resolved structure. Secondly, its implied size calculated as 
a geometric thin-ring radius for NGC4151 as well as those for 
the other type 1 s are found to roughly match independent radius 
measurements from near-IR reverberations, which are thought to 
be probing the dust sublimation radius. 

The near-IR reverberation radius Rj^ is the light-crossing 
distance over the time lag between the variabilities in the 
UV/optical and near-IR (2.2 fim). The former emission is 
thought to originate from the central accretion disk (AD) and the 
latter from the hot dust grains in the sublimation region heated by 
the AD (e.g.lG lass 2004; M inezaki et al.ll2004t ISuganuma et al] 



1.3. Type 1 study in the mid-IR 



l2006t [Koshida et aL, ,2009.) . The radius R^^^ has been shown 
to be p roportional to L'^^ where L is the AGN luminosity 
( ISug anuma et al. 2006; our specific definition of L is given in 
Sect l3.ll ). thus, the interferometric ring radii at K-band described 
above also roughly scale with L^^^. This is what we w ould expect 
for the dust sublimation radii (e.g. lBarvaini i[l983). In addtion, 
the measured ring radii have been shown to be almost unaffecte d 
by the near-IR AD component (Kishimo to et alJl2007l l2009ah . 
Therefore, the case for the near-IR visibiUties probing the dust 
sublimation region is strongly made. 

Furthermore. lKishimoto et all (l2009aLl20Tll) pointed out that 
the ring radii derived from the observed visibilities are either 
roughly equal to or slightly larger than the L'^^-fit to the rever- 
beration radii as a function of L. Since reverberation measure- 
ments general ly place weight on the small responding radii (e.g. 
iKoratkar & Gaskelllll99ll) . Rt^^ is expected to represent the ra- 
dius very close to the inner boundary of the dust distribution. On 
the other hand, the interferometric ring radii rather correspond 
to brightness-weighted effective radii. In this sense, the L'^^-fit 
reverberation radius gives a good normalization for the radial 
extent, taking out the simple 1}^^ scaling. This normalization is 
used extensively in this paper 



In the mid-IR, iBeckert et all ( l2008h presented the first long- 
baseline mid-IR interferometry for a type 1 AGN, NGC3783. 
Kishimoto et al. (2009b) published two-ba seline data f or the 
same galaxy, revising the results shown in iBeckert et al. l (E008l) 
with dedicated software (see more details below) and also adding 
one more data set taken at a different baseline. They proposed a 
generic way of implementing a uniform comparison of the inter- 
ferometric data from different targets with different luminosities 
and at different distances, by using the inner radius Ri^ described 
above for normalizing the spatial scale probed by the inteferome- 
ter. Mainly using the data of NGC3783, but supplementing them 
with those from other type Is and the type 2 AGN NGC1068, 
they argued that the radial brightness distribution in the mid-IR 
of type 1 AGNs seems to follow a relatively shallow power-law. 
This was a first attempt to directly map the radial distribution in 

type 1 AGN tori. 

Subsequently, iBurtscher et al.l (l2009l) also presented two- 
baseline MIDI data for another type 1 AGN, NGC4151, and 
measured the mid-IR emitting region size using a model of a 
Gaussian plus point-source, which gives a visibility curve rela- 
tively similar t o th at of a power-law at low spatial frequencies. 
iTristram et al.l (l2009i) presented the results of a MIDI snap-shot 
survey for AGNs of both types 1 and 2, and discussed the mid-IR 
emission size as a function of luminosit y, which has recently be 
further evaluated with additional data bv lTristram & SchartmannI 
(201J1). 

Here we present a first systematic investigation of type 1 
AGNs with the mid-IR interferometer, expanding the study of 
Kishimoto et al. (2009b). Combining our mid-IR data with the 
near-IR interferometry and photometry, we aim to comprehen- 
sively map the radial structure of the thermally emitting dusty 
region in the mid- and near-IR. We aim to derive physical con- 
straints as directly as possible from the observations, with the 
ultimate goal of finding whether the radial structure is systemat- 
ically dependent on the properties of the central engine. 



2. Observations and data reductions 

2.1. MIDI 

We observed six type 1 AGNs listed in Table [1] on 8-12 May 
2009 (UT) using MIDI, wh ich is a two-beam c ombiner working 
in the mid-IR wavelengths dLeinert et al.ll2003h . The observation 
log for each fringe track is listed in Table |2] We also collected a 
limited number of data sets for two of the six targets on 15 Mar 
and 5 Aug 2009 (Table|2]i. The adaptive optics system MACAO 
at each of the 8.2m Unit Telescopes was locked on the nucleus of 
each target AGN. A single baseline was allocated for each night, 
and we used four different baselines with their projected length 
spanning from 27 m to 129 m. For one of the targets, NGC3783, 
we also used in our ana lysis MIDI data taken in 2005 published 
bv .Beckertet all ( l2008h and Kishimoto et al. (2009b). 

T he data wer e reduced using the software EWS (version 
1.7.1: iJaffd [2004 . but with modifications using our own IDL 
codes. The old NGC3783 data were also re-reduced with the 
same procedure. The details of the process are described in the 
AppendixlAl Briefly, we firstly implemented an additional back- 
ground subtraction for the photometry frames using adjacent sky 
strips in the 2D spectrum. Secondly, to reduce errors in group 
delay determinations, we smoothed delay tracks over typically 
~20-40 frames. To avoid possible positive biases in the corre- 
lated flux, we also averaged ~20-40 frames to determine the 
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Table 1. Properties of our targets and their point-source flux from our UKIRT images or 2MASS images. 



name 




scale^ 




date 






lux (mJy) 






Av 


Sr^ fit^ 




yLi,(55uOA)' 




corr. 


(pc/mas) 


(mag) 


(UT) 


Z 


Y 


J 


H 


K 


(mag) 


(pc) 


(mas) 


(ergs-')* 


NGC4151 


0.00414 


0.0855 


0.028 


2009-06-17 


39.7 


42.4 


59.4 


96.3 


173. 




0.033 


0.383 


(8.3 ± 1.3) ■ lO"*^ 


NGC3783S 


0.0108 


0.221 


0.119 


1999-04-11 






19.1 


33.5 


62.5 




0.061 


0.274 


(2.8 ± 0.4) ■ 10"' 


ES0323-G77'; 


0.0159 


0.324 


0.101 


1999-05-03 






28.6 


63.4 


107. 


1.2±0.4 


0.10 


0.321 


(8.4 ± 1.9)- lO"*' 


H0557-385 


0.0342 


0.681 


0.043 


2009-12-13 


4.03 


4.43 


7.21 


15.3 


36.8 


1.5±0.5 


0.12 


0.181 


(1.2±0.3)- lO""* 




0.0579 


1.12 


0.182 


2000-01-05 






35.3 


59.1 


116. 




0.47 


0.419 


(1.7±0.2)- lO''^ 


IRAS 13349' 


0.108 


1.98 


0.012 


2009-06-17 


7.75 


8.71 


12.6 


26.3 


62.3 


0.93±0.31 


0.49 


0.248 


(1.9±0.4)- lO"*^ 



" CMB coiTected value from NED . * //p = 70 km s ' M pc n„, = 0.3, and = 0.7. ' Galactic reddening from Schlegel et al. 1 1998). Flux measuremen t 
uncertainty is ~5 % (see Sect.3 in iKishimoto et al]|2009al) . ' The K-band reverberation radius Rj^ from UV luminosity using the fit bv Suganuma et al] j2006h . 
which is defined as Ri„ in this paper (see Eq[T). ^ Optical luminosity at rest wavelength 5500A. From 2MASS images. IRAS09 149-6206. ' IRAS13349+2438. 



phase offsets. The system visibihty was obtained from the ob- 
servations of visibihty cahbrators taken right after or before each 
target observation with a similar airmass. These data were also 
reduced with similar smoothings above to calibrate out the effect 
of the time averaging. These calibrators are also the photometric 
standards found in the list by C ohen et al.l (I1999D and/or by R. 
van Boekel (priv. communication). We obtained the correlated 
flux, which was corrected for the system visibility, and total flux 
separately. The errors for the total flux and correlated flux were 
firstly estimated from the fluctuation of the measurements over 
time. The total flux spectrum for each object was then derived as 
the weighted mean of many measurements over diff'erent nights 
(for NGC3783, the spectrum was derived separately for 2005 
and 2009 data; see below). For the correlated flux measurements 
at adjacent uv points, we also took weighted means. Finally, for 
the correlated flux, we assigned a systematic uncertainty as a 5% 
of the total flux based on the analysis presented in Appendix lAl 

For NGC3783, the mid-IR total flux apparently increased 
by ~20% from 2005 to 2009 observations. To facilitate a di- 
rect comparison, we scaled the total flux and correlated flux 
from the 2005 observations by the same factor This means 
that we assumed no visibility change in the mid-IR over this 
period. This is conceivable based on a comparison of near- 
IR interferometry data for NGC4151 that fou nd no significant 
visibility change over a one-year timescale (iPott et al.l 1201 Ot 
[Kishimoto et al. 2011). 

To test the accuracy of our MIDI reduction software specifi- 
cally for targets at sub-Jy flux levels, we observed stars with flux 
~0.3-0.8 Jy at 12 fim during our May 2009 run, and also in Mar 
2009. Details are described in Appendix lAl One significant is- 
sue is that, if the frames are not averaged over a sufficiently long 
time-interval, the correlated flux is affected by a positive bias 
apparently in the wavelength region for which there are lower 
(intrinsic) counts. Since MIDI count spectra often have a peak at 
shorter wavelengths, the resulting biased spectrum is often red- 
der than it should be. In addition, coherence loss toward shorter 
wavelengths might not be fully compensated for by calibrator 
frames, leading also to redder correlated-flux spectra. With the 
data for the sub-Jy unresolved stars, we demonstrate that the 
dedicated codes can correctly estimate the correlated flux, and 
quantify the possible systematic uncertainty. 

In Table |2] the quality of each fringe track is indicated by 
the maximum signal-to-noise ratio (S/N) per smoothed frame 
(see Sect lA.3b . The analysis of the sub-Jy unresolved stars de- 
scribed above currently covers this S/N down to ~2.5, thus in 
this paper we decided not to use data with a S/N lower than 2.5. 
Accordingly, five fringe tracks in Table |2] were excluded from 
our analysis. 



2.2. The near-IR data 

We combined the mid-IR MIDI data above with data in the 
near-IR to comprehensively map the radial structure of the ther- 
mally emitting dusty region. For the photometry in the near- 
IR, we have contemporaneous imaging data at five broad-band 
wavelengths from 0.9 to 2.2 fim for NGC4151, IRAS13349, 
and H0557-385, taken with the Wide Field Camera (WFCAM) 
on the United Kingdom Infrared Telecope (UKIRT ). The data 
for the former two objects were published in Kish imoto et al.l 
( l2009al) . but, together with the new data for H0557-385, we 
present here the data for all three AGNs in Table [1] The wide 
field-of-view of WFCAM provides simultaneous point-spread- 
function (PSF) measurements from stars in the same field. 
We used them to implement accurate PSF + host galaxy fits, 
and measured the poin t-source-only SEDs in the near-IR (see 
IKishimoto et al ] |2009al for more details). For the three other tar- 
gets in the sample, we measured the point-source flux by apply- 
ing the same procedure to 2MASS images. 

We have the near-IR interferometry data taken with the Keck 
interferometer (KI) at 2.2 jum for two of the six MIDI targets, 
NGC 4151 and IRAS 13349. T he KI data were taken on 15 May 
2009 ( IKishimoto et al.ll2009a \ less than a week from our MIDI 
observations of these targets. For the general comparison over 
the sample between the mid- and near-IR, we also use the sample 
of in tot al eight type 1 AGNs, taken w ith the KI in the near-IR at 
/./ /im (IKishimoto et alJl2009aLl201 ih . 



2.3. VISIR 

To obtain accurate photometry and check the accuracy of the to- 
tal flux spectra from MIDI, we took mid-IR images of each ob- 
ject in two passbands, namely at 8.6 yum and 11.9 fim with VISIR 
on UT3. These images were mostly taken within about four 
weeks of the MIDI observations, to minimize the effect of vari- 
ability. The observation log and the results are show n in Table |3l 
For ES0323-G77, a VISIR spectrum was taken bv lHonig etall 
(j20^10ft within several days of our MIDI observations. 

All the VISIR images were taken in 'perpendicular' mode 
with each dataset consisting of four image tiles. Using the data 
from the ESO pipeline (version 3.7.2), we implemented addi- 
tional subtraction of the background determined at annuli of 
radii ~ 1.5 - 3.0" from each source on each of the four tiles. 
After checking the stability of the flux calibration as a function 
of aperture radii, we determined the flux conversion factor at a 
radius ~l" where its uncertainty was estimated from the scatter 
in the results from the four tiles. 
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Table 3. The results of our VISIR imaging observations. 



date & time 


telescope 


Bp 


PA 


vis./flux 


max 


(UT) 






(m) 


(°) 


calibrator 


S/N° 


NGC4151 














2009-05-09 


01:20 


UT1-UT3 


61.5 


55.9 


HD94264 


4.1 


2009-05-09 


01:32 


UT1-UT3 


64.1 


55.7 


HD94264 


3.3 


2009-05-09 


02:37 


UT1-UT3 


77.3 


53.0 


HD98262 


3.3 


2009-05-09 


02:49 


UT1-UT3 


79.3 


52.3 


HD98262 


3.4 


2009-05-10 


00:44 


UT1-UT2 


26J 


45.1 


HD94264 


5.1 


2009-05-10 


00:55 


UT1-UT2 


28.1 


46.1 


HD94264 


4.9 


2009-05-11 


00:33 


UT1-UT4 


102.0 


84.8 


HD 113996 


4.1 


2009-05-11 


00:44 


UT1-UT4 


105.1 


83.1 


HD 113996 


4.3 


2009-05-11 


01:19 


UT1-UT4 


113.4 


78.3 


HD94264 


4.2 


2009-05-11 


01:31 


UT1-UT4 


115.8 


76.7 


HD94264 


4.0 


2009-05-12 


01:59 


UT1-UT4 


121.3 


72.4 


HD94264 


3.6 


2009-05-12 


02:28 


UT1-UT4 


125.2 


68.3 


HD98262 


3.3 


NGC3783 














2005-05-28 


03:39 


UT1-UT2 


44.4 


45.1 


HD 100407 


3.1 


2005-05-28 


04:00 


UT1-UT2 


42.6 


45.9 


HD 100407 


5.6 


2005-05-31 


03:01 


UT2-UT4 


68.7 


114.8 


HD 100407 


3.2 


2005-05-31 


03:21 


UT2-UT4 


64.9 


119.6 


HD 100407 


3.7 


2009-05-11 


03:31 


UT1-UT4 


108.3 


80.4 


HD112213 


3.7 


2009-05-11 


23:19 


UT1-UT4 


128.4 


45.0 


HD101666 


4.5 


2009-05-11 


23:33 


UT1-UT4 


129.0 


47.3 


HD101666 


4.2 


ES0323-G77 














2009-05-08 


01:01 


UT3-UT4 


56.8 


94.4 


HDl 12213 


2.1* 


2009-05-08 


01:06 


UT3-UT4 


57.3 


95.3 


HDl 12213 


2.5* 


2009-05-10 


03:06 


UT1-UT2 


54.1 


30.5 


HD112213 


3.8 


2009-05-10 


03:18 


UT1-UT2 


53.8 


31.8 


HD112213 


3.4 


2009-05-10 


04:09 


UT1-UT2 


51.8 


37.1 


HDl 12213 


3.8 


2009-05-10 


04:21 


UT1-UT2 


51.3 


38.1 


HD112213 


4.1 


2009-05-11 


04:07 


UT1-UT4 


118.4 


75.0 


HD112213 


4.2 


2009-05-12 


00:05 


UT1-UT4 


126.9 


36.9 


HDl 12213 


3.9 


2009-05-12 


00:16 


UT1-UT4 


127.6 


39.1 


HD112213 


3.6 


H0557-385 














2009-05-09 


23:33 


UT1-UT2 


42.1 


46.4 


HD40808 


3.6 


2009-05-09 


23:45 


UT1-UT2 


41.0 


46.8 


HD40808 


4.0 


2009-08-02 


09:14 


UT3-UT4 


32.9 


60.3 


HD9362 


3.8 


2009-08-03 


09:02 


UT1-UT2 


56.3 


-12.9 


HD39425 


3.7 


2009-08-05 


09:48 


UT1-UT4 


120.0 


17.4 


HD40808 


3.5 


IRAS09149-6206 












2009-03-15 


00:29 


UTl-tJT2 


49.9 


13.0 


HD80007 


3.9 


2009-03-15 


00:39 


tJTl-tJT2 


49.7 


14.7 


HD80007 


3.8 


2009-05-08 


00:25 


UT3-UT4 


62.4 


126.4 


HD89682 


2.3* 


2009-05-08 


01:45 


UT3-UT4 


62.4 


144.2 


HD89682 


2.5* 


2009-05-08 


02:01 


UT3-UT4 


62.3 


148.0 


HD89682 


2.0* 


2009-05-09 


00:18 


UT1-UT3 


80.2 


53.4 


HD89682 


3.6 


2009-05-09 


00:30 


UT1-UT3 


78.9 


55.4 


HD89682 


4.2 


2009-05-10 


23:38 


UT1-UT4 


119.2 


77.4 


HD89682 


4.2 


2009-05-10 


23:49 


UT1-UT4 


118.0 


79.8 


HD89682 


4.6 


IRAS13349-(-2438 












2009-05-10 


01:51 


UT1-UT2 


34.7 


28.6 


HDl 13996 


3.2 


2009-05-10 


02:30 


UT1-UT2 


38.2 


33.6 


HDl 13996 


3.2 


2009-05-10 


04:40 


UT1-UT2 


49.3 


39.3 


HDl 12213 


4.6 


2009-05-11 


02:40 


UT1-UT4 


116.2 


71.1 


HDl 13996 


4.5 


2009-05-11 


02:51 


UT1-UT4 


118.7 


70.3 


HDl 13996 


4.9 


2009-05-12 


01:06 


UT1-UT4 


89.9 


76.2 


HDl 13996 


3.5 


2009-05-12 


01:18 


UT1-UT4 


93.9 


75.6 


HDl 13996 


3.3 


2009-05-12 


04:11 


UT1-UT4 


129.3 


63.7 


HDl 27093 


5.1 



" Maximum signal-to-noise ratio (SNR) per smoothed frame for the correlated 
flux integraged over the whole N-band (see Appendix for details). 
* Excluded from the analyses in this paper due to low SNR. 



2.4. Additional piiotometric data in tine mid-IR and near-IR 

When fitting the models described in the next section to the 
SEDs, we only used our MIDI total flux spectra and WFCAM 
(or 2MASS) point-source fluxes described above. However, we 
gathered in addition IRS spectra and IRAC imaging data from 
the Spitzer Heritage Archive, listed in Table |4] for all the tar- 
gets. For NGC3783, we also have VLT/ISAAC data taken on 2 
Jan 201 1 (Honig et al. in prep; Table|4]i. These extra data points 
from Spitzer and ISAAC fall approximately on our model SEDs 



name 


flux 


aperture" 


UT date 


calibrator 




(mJy) 


(arcsec) 






PAHl at A 8.59 fim, AA 0.42 






NGC4151 


860±17 


1.8 


2009-06-07 


HD120933 


NGC3783 


423±24 


1.5 


2009-06-07 


HD 145897 


H0557-385 


330±18 


1.5 


2009-09-07 


HD41047 


IRAS09149 


346±10 


1.5 


2009-06-07 


HD99167 


IRAS 13349 


389±10 


1.5 


2009-06-07 


HD145897 


PAH2ref2 at A 11.88 yum, 0.37 yum 






NGC4151 


1432±53 


2.7 


2009-06-07 


HD120933 


NGC3783 


760+17 


1.5 


2009-06-07 


HD145897 


H0557-385 


404+29 


1.5 


2009-09-07 


HD41047 


IRAS09149 


474±31 


1.5 


2009-06-07 


HD99167 


IRAS 13349 


530±26 


1.5 


2009-06-08 


HD145897 


" Aperture diameter for target and calibrator for flux determination. 



presented below, with possibly some additional flux from the 
host galaxy, thus are viewed as approximate upper limits, and 
not used in the spectral fitting. We summarize the usage of our 
data in Table |5] 



2.5. Reddening corrections 

We corrected t h e S ED for Galactic reddening using values from 
ISchlegel et aT! (1998) hsted in Table □ We used the redden- 
ing ^curves_ofJ^Mdellietal] (11989) . supplemented with that of 
Chi ar & TielensI (120061) at long wavelengths. We also corrected 
the SED for the possible reddening within the host galaxy for 
IRAS 13349, ES0323-G77, and H0557-385, using the Ay val- 
ues listed in Table [T] For H0557-385, in terms of the optical 
continuum slope, Balmer decrement, and host galaxy inclination 
(Winkler 1992; Fairall et al. 1982), the reddening is likely to be 
between those of ES0323-G77 and IC4329A, which is another 
Seyfert 1 galaxy known to be strongly affected by host redden- 
ing. For the latter two objects, the estimates of Ay from the dif- 
ferent methods exist including that derived from the gradients of 
flux va riation ( Winkler et al. 1992), and are ~1 and ~2, respec- 
tively (IWinkler et al.lll992D . Without any additional solid infor- 
mation about the reddening value for H0557-385, we roughly 
infer that Ay = 1.5 +0.5. 



3. Results 

3.1. The whole data sets and prerequisites 

Figures [1] and |2] present all the data for each of the six targets. 
The total and correlated flux spectra are shown in the top panel, 
the visibility as a function of spatial frequency in the middle, 
and the uv coverage in the bottom. In the top panel, all the addi- 
tional flux data described in the previous section are also plotted. 
The MIDI data were binned by seven pixels (AA ~ 0.4jum) with 
the error in each bin taken as the median of the errors over the 
binned spectral channels (since eiTors over adjacent channels are 
coiTelated). 

For the visibilities, the spatial frequency is shown in units 
of cycles per Rin. As described in Sect. 11.21 as the most empir- 
ical measure of the inner radius of the dust distribution, we 
ado pt the L^^^ fit to the K -band reverberation radius Rr^, given 
by Suganum a et al.l (l2006h . Here L is the UV luminosity of the 
central engine, which is defined as a scaled optical luminosity, 
L = 6vLv(5500A) (iKishimoto et al. 2007), assuming a generic 
SED over the UV/optical wavelengths. Here v designates the 
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Table 4. Spizer data from Spitzer Heritage Archive and ISAACA'^LT data. 



name 


IRS 






IRAC flux (mJy) 










date (U 1) 


date (UT) 


3.6 yum 


4.5 fim 


5.7 


7.9 yum 




NGC4151 


2004-01-08 


2004-12-17 


264.9±1.0 


363.8+1.2 


554.0+3.2 


924.7+2.3 




NGC3783 


2007-06-27 


2007-08-07 


118.0+0.7 


158.5+0.8 


228.1+2.3 


359.6+1.5 




ES0323-G77 


2006-08-05 


2009-03-14 


185.6±0.8 


215.4+1.0 


288.3+2.4 


411.2+1.6 




H0557-385 


2007-10-06 


2008-10-31 


I40.7±0.7 


195.3+0.9 


268.8+2.4 


344.0+1.5 




IRAS09149 


2008-07-02 


2008-05-13 


214.4±0.9 


261.9+1.1 


330.9+2.6 


391.8+1.6 




IRAS 13349 


2005-06-07 


2009-02-02 


199.1±4.2 


242.6+1.0 


300.1+4.2 


431.4+1.7 




name 




ISAAC flux (mJy) 












date (UT) 


2.25 yum 


3.80 yum 


4.66 yum 








NGC3783 


2011-01-02 


66.21±0.61 107.00±0.49 152.83+4.65 








Table 5. Summary of the data used for each target. 










name 




data used for modeling 




additional photometric data only for 


polar- 






near-IR 




mid-IR 




confirmation or approximate upper limits 


axis" 




1-2.2 ycim flux interferometry 


total flux interferometry 


3-15 yum mid-IR 


PAD 


NGC4151 


WFCAM/UKIRT Keck 


MIDI spectrum 


MIDI 


IRAC / IRS VISIR images 


91" 


NGC3783 


2MASS 




MIDI spectrum 


MIDI 


ISAAC / IRAC / IRS VISIR images 


136" 


ES0323-G77 


2MASS 




MIDI spectrum 


MIDI 


IRAC / IRS VISIR spectrum 


174'' 


H0557-385 


WFCAM/UKIRT 


MIDI spectrum 


MIDI 


IRAC / IRS VISIR images 


33" 


IRAS09149 


2MASS 




MIDI spectrum 


MIDI 


IRAC / IRS VISIR images 




IRAS 13349 


WFCAMAJKIRT Keck 


MIDI spectrum 


MIDI 


IRAC / IRS VISIR images 


35-^ 



"The PA of the polar axis of the system inferred from optical polarization and radio axis. 

iMarteL Cl99&):,Mundell et al.. l2003j. '.Smith et al.. (,2002). '^.Schmid et al.. C2003). '.Brindle et alJ llI99(]h. IWills et all lll992h. 



frequency in the rest frame, and Ly is the isot ropic luminosity 
per fre quency. T he formula for the L'^^ fit by ISuganuma et alj 
(l2006l) is given in lKishimoto et all (l2007h as 



Ri, = = 0.47 ■ 



6vLv(5500A) \ 
10"^*' erg /sec / 



1/2 



pc. 



(1) 



If the dust temperature structure is primarily determined by 
the illumination from the central source, the radial temperature 
structure scales with the sublimation radius In this case, by 
normalizing a length scale with j^jn, we can place various inter- 
ferometric data on a uniform scale (iKishimoto et al.ll2009b) . A 
given configuration of an interferometer probes a certain spa- 
tial frequency, or a certain spatial scale represented by a spa- 
tial wavelength, which is the reciprocal of the spatial frequency. 
Here we normalize this spatial wavelength with Ri^ (see the up- 
per axis of the middle panel of FigslT]and|2|. This corresponds 
to writing the spatial frequency in units of cycles per Ri^. With 
this normalization, we can eliminate the simple luminosity scal- 
ing originating from R^^ being proportional to L''^, as well as 
a simple distance scaling. Different objects with different lumi- 
nosities and at various distances have different angular sizes for 
Rin, but the data can be uniformly viewed when the probed spa- 
tial scale is normalized with Ri^. We can further directly deter- 
mine whether there is a luminosity dependence in the structure 
beyond the simple L'^^ scaling. 

Following the same idea, the total and correlated fluxes can 
be normalized using R^^. They are shown in the top panel of 
Figs[T]and|2]as {vLy I An) I r1^^ and have the same dimension as the 
surface brightness in the rest frame of each target. This is equal 
to Vobs/y„b, •(! where Vobs, /y„b,, and ^i,, are the frequency, 

observed flux, and angular size of in the observed frame, 
respectively. The correlated flux in the top panel is color-coded 
with the spatial wavelength, while each visibility data point in 
the middle panel is color-coded with the (radiation) wavelength 



in the rest frame, as indicated by the color bars. The total flux, 
correlated flux, and visibility spectra are shown in Appendix IB] 
using more conventional units to facilitate direct comparisons 
with the data in the literature. 



3.2. MIDI total flux spectra 

With MIDI, total fluxes are measured separately and indepen- 
dently of correlated flux measurements. We derived a total flux 
spectrum for each object as a weighted mean of many measure- 
ments over different nights. The fluxes from our VISIR photom- 
etry (or VISIR spectrum in the case of ES0323-G77) were found 
to be either in excellent agreement (within a few percent) with 
or slightly larger (~10-20%; see Appendix IbTi than the mean 
MIDI total flux spectrum for each object. Since the latter is from 
the beams corrected with MACAO and from a smaller extrac- 
tion window (adaptive mask in EWS), it is conceivable that it 
has less contribution from the extended emission, e.g. from the 
host galaxy. In the case of NGC4151, which seems to have quite 
an extended mid-IR emission (see the low visibilities in Fig[T^), 
some of the extended emission in the VISIR flux is probably 
excluded in the MIDI flux. All these comparisons with VISIR 
fluxes confirm the accuracy of the mean MIDI total flux, and its 
error might be slightly over-estimated in some cases. 

In the middle panel for visibilities in FigslT]and|2](and also in 
FigHdiscussed in the next section), we show the visibility uncer- 
tainties originating from each correlated flux measurement with- 
out the contribution of total flux errors, but show the total flux 
errors separately as data points of visibility unity at the lowest 
spatial frequency. In this way, we can clearly show the relative 
uncertainties between each correlated flux measurement and the 
total flux measurement. In addition, these total flux data points 
are shown at the spatial frequency with the baseline length equal 
to the single telescope diameter. As long as the source remains 
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color-coded on spotiol wavelength (Ri„) color-coded on spotiol wovelength {R,„) color-coded on spotiol wavelength (R^) 




n n r . . . r. , ,.1 n n r . . ? n n r . . ^^^^ 

0.001 0.010 0.100 0.001 0.010 0.100 0.001 0.010 0.100 

Spatial frequency (cycles per R|„) Spatiol frequency (cycles per R;„) Spatial frequency (cycles per R|„) 



coloi — coded on rest wavelength O^m) color-coded on rest wavelength (;im) color— coded on rest wavelength {f^r^) 




U (m) U (m) U (m) 

(a) NGC4151 (b) NGC3783 (c) ES0323-G77 

Fig.l. The observed data for (a) NGC4151, (b) NGC3783, and (c) ES0323-G77 (see SectHI, and the best fits using a tempera- 
ture/density gradient model plus an inner 1400K ring (see Sect l4.4l and Table[8l). Top: Spectral energy distribution in {vLy/4jT)/R^^. 
The data plotted in black are the total flux measurements from MIDI and VISIR in the mid-IR and from WFCAM or 2MASS in the 
near-IR. Additional total flux spectra and photometry from Spitzer (Table|4]l are plotted in gray. For NGC3783, ISAAC photometry 
data (Table m are shown as gray squares. Correlated fluxes from the mid-IR and near-IR interferometry are plotted with colors 
coded on spatial wavelengths in units of R\„ as shown in the top color bar. The model curves for the total flux and the correlated 
fluxes at baselines corresponding to each MIDI observation are drawn in the same color coding. Middle: Visibility amplitude plotted 
against spatial frequency in units of cycles per Ri„. The color of the data points corresponds to the observing wavelength in the rest 
frame of the target. This color-coding is shown as a color bar between the middle and bottom panels. The iso-wavelength model 
visibility curves are shown for 13, 8.5, and 2.2/im, as well as the curves for each baseline configuration (which are 'orthogonal' 
to the iso-wavelengh curves). Dotted curves are the visibilities only for the power-law component, with the difference from solid 
curves indicating the effect of the near-IR ring. For NGC3783, solid and dashed model curves are for equatorial and polar axis 
PAs, respectively (see Sect l4.3l l. Bottom: The uv coverage of the mid-IR and near-IR interferometry with the same color-coding on 
Observing wavelengths. Solid and dashed lines indicate the PA of the projected equatorial and polar axes, respectively, inferred from 
optical polarization and radio data (see Table|5]l. 
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(a) H0557-385 (b) IRAS09149 

Fig. 2. The same as FigH] but for (a) H0557-385 (b) IRAS09149 (c) IRAS 13349. 



150 100 50 -50 -100 -150 
U (m) 

(c) IRAS 13349 



unresolved by a single telescope, the total flux gives the corre- 
lated flux at this spatial frequency, and the visibility is unity at 
this spatial frequency. Thus, the data points illustrate the con- 
straint in such a case, which is expected to hold at least approx- 
imately for our sample of type 1 objects. We note however that 
these are only for the data presentation. For the model fitting in 
Sect l3.5l and [H we include the total flux error in the visibility 
error budget in a normal way. 



3.3. Mapping radial structure 

Each panel of Figs[T]and|2]shares a common range of axes over 
different objects so that we can directly compare them. We can 
also place together the data for all the targets on a single plane. 
Fig|3]is such a figure, showing the observed mid-IR visibilities 
of all the targets, each averaged from 10 to 12 jim in the rest 
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Fig. 3. Observed mid-IR visibilities of all the targets, each averaged over 10 to 12 fim in the rest frame of each object, shown as 
a function of spatial frequency in units of cycles per Ri„ (log-scale). The data are color-coded in terms of the UV luminosity of 
the central engine for the corresponding object, which is defined to be 6 vLv(5500A) dKishimoto et alj|2007h . Each data point is 
plotted with different symbols for different objects, and is labeled with the PA / length of the baseline. It is quite clear that the radial 
structure in units of R^n changes from being compact in higher luminosity objects to extended in lower luminosity objects. Dotted 
curves are visibility functions for power-law brightness distribution with the normalized half-light radius Ri/2/Rm-3, 10, and 30., 
which roughly depict the brightness distributions from higher to lower luminosity objects. 



frame of each object. In Appendix ICl we also show the same 
figure with the spatial frequency axis on a linear scale. 

Our mid-IR observations cover the range of several tens to 
hundreds of Ri„ in spatial wavelengths. Over this range, the 
observed visibilities are relatively concentrated in one locus, 
from 0.3 to 0.9. If the visibilities of different objects followed 
the same curve, this would mean that the size scales only with 
L^^^ a nd objects share the same radial structure as inferred by 
iKishimoto et al.l (l2009bl) . However, this is not the case. 

In Fig. |3] each object is color-coded with its UV luminosity 
as defined in Sect l3.ll It is quite clear from the figure that within 
the small sample, which nevertheless spans over ~2.5 orders of 
magnitudes in luminosity, the radial structure does depend on 
luminosity, i.e. the structure in units of Ri^ looks more compact 
in higher luminosity objects. This implies that the increase in the 
mid-IR physical size (e.g. in pc) with luminosity is more gradual 
than L'^^. We discuss this luminosity dependence in Sect l4.1l 

While the results above are essentially model-independent, 
we attempt in the following to describe the radial surface- 
brightness distribution more specifically with simple models. 
The dotted curves in Fig. [3] are visibility functions for power- 



law brightness distributions, which we argue below approxi- 
mately depict the observed brightness distributions (see more in 
SectlTSll. 

3.4. Single Gaussian/ring versus power-law 

We first consider the results for ES0323-G77, where we have 
the data with the least ambiguity. The total flux spectrum pro- 
vides with an excellent match to the simultaneous VISIR spec- 
trum (see Sect j2.3l Fig jB.ll i and we have the two datasets at al- 
most exactly the same PA of the baselines (~38°) with good S/N. 
The visibility data are shown in Fig|4]as a function of spatial fre- 
quency. 

As in the middle panel of Figs. [T]and|2l the visibilities in 
Fig|4]are presented in different colors for different wavelengths 
as indicated by the color bar. For a given baseline, shorter wave- 
lengths probe the source at higher spatial resolutions, and if the 
overall source size is the same when seen at different wave- 
lengths, then the visibility is expected to become lower, i.e. more 
resolved, toward shorter wavelengths. However, the observed 
visibility spectra are roughly flat, or increase toward shorter 
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Fig. 4. Top: Observed visibilities for ES0323-G77 as a function 
of spatial frequency along the position angle ~38°. Visibility 
curves for Gaussian distributions with a HWHM//?in of 20 
(dotted) and 12 (dashed), and a power-law distribution with 
RililR'm - 7 (dash-dot) are over-plotted. Bottom: Gaussian 
HWHM (plus signs) and R\i2 (filled circles) are plotted against 
the same spatial frequency. The red to green colors indi- 
cate the observing wavelengths in the rest frame from 13 to 
8.5 yum. For a given wavelength, while the Gaussian HWHM 
is baseline-dependent, the half-light radius is approximately 
baseline-independent, thus seems adequate for representing the 
whole structure. 

wavelengths, which means that the overall size of the source 
is smaller at shorter wavelengths. This is also generally seen 
in other objects, except at some objects' shortest wavelengths 
where visibilities might be slightly underestimated owing to a 
coherence loss (see Appendix lA.4b . 

On the other hand, when we consider the visibilities at each 
wavelength as a function of baselines (or spatial frequencies), 
we can easily see that a simple Gaussian or ring (which we treat 
here almost equally, since they give almost identical visibility 
curves at low spatial frequencies well before the first null; see 
Sect l4.2l ) does not fit them well, as shown by the dotted and 
dashed Gaussian curves in the top panel of Fig|4] If we derive 
a Gaussian HWHM or ring radius for each baseline measure- 
ment at a fixed wavelength, we obtain a smaller size for longer 
baselines as we show in the bottom panel of FiglD We appear 
to observe here a brightness distribution where the correspond- 
ing Gaussian/ring sizes become progressively smaller at longer 
baselines. This probably means that there is a continuous distri- 
bution of brightness spread over a wide range of continuous spa- 
tial scales. This can be understood if the brightness distribution 
resembles a power-law distribution at a given wavelength, with 
its index generally becomi ng steeper at shorter w avelengths. 
This is what we advocated in lKishimoto et al ] (l2009b ) (see their 
Fig.l). This power law is at least partially expected physically, 
since the temperature of the directly illuminated dust grains is 
expected to follow a power-law with the radius from the il- 
lumination sou rce, simply based on thermal equilibrium (e.g. 
lBarvainisl[T987h . 



0.01 0.10 

Spatial frequency (cycles per Rj/j) 

Fig. 5. Top: Comparison between visibility functions for 
Gaussian, shown as a dashed curve, and power-law distributions 
where red, green, and purple curves correspond to the cases of 
R\I2 - 6, 2, 1.3 Rin, respectively. For comparison, the visibility 
function of a thin ring with radius being equal to the Gaussian 
HWHM (i.e. Rxji) is shown as a dotted curve. Bottom: The ratio 
of Gaussian HWHM to /?i/2 is plotted against the same spatial 
frequency in each case. The vertical dot-dash line corresponds 
to the spatial wavelength of 4.5 Riji- 



3.5. Half-light radius for power-law brightness distribution 

For given inner and outer radii, a power-law index specifies a 
power-law brightness distribution, and a convenient way of char- 
acterizing the corresponding size is to use the half-light radius 
R\/2, the radius within which half of the total integrated light 
at a given wavelength is contained. This can be directly com- 
pared with the HWHM of a Gaussian brightness distribution, 
since R1/2 is equal to the HWHM in a Gaussian case. When the 
brightness distribution is adequately described by a power law, 
the values of 7? 1/2 derived from visibilities at different baselines 
should all be the same for a given wavelength. The bottom panel 
of Fig m presents the deduced R1/2 for each visibility observed. 
We can see that the deduced R1/2 is quite independent of the 
baseline lengths, in contrast to the Gaussian HWHM. This re- 
assures that the power-law description in this case is quite ad- 
equate. We can also see that the size becomes smaller toward 
shorter wavelengths, as discussed above. 

The disadvantage of adopting a power-law model and its cor- 
responding Ri/2 is that we need to define both the inner and 
outer radii, which are set here to be IRi^ and 1007?in, respec- 
tively. However, we can also define another, very similar, more 
model-independent radius, the half-visibility radius 7?vo.5 as fol- 
lows. A visibility at a given low spatial frequency (within the 
first lobe of the visibility function before the first null) roughly 
corresponds to the fraction of the flux that is contained within 
the spatial scale being probed and remains unresolved by the 
interferometer's configuration. Here the configuration is set by 
the observing wavelength A and baseline length B. In this sense, 
Ri/2 roughly corresponds to the spatial scale probed by the in- 
terferometer for which the visibility becomes 0.5. We can re- 
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Fig. 6. Normalized half-light radii R\/2/Rm plotted against rest 
wavelengths. Thin-ring radii a t 2.2^m, normali zed by Rin, from 
the KI data for two objects dKishimoto et alj l2009al) are also 
shown. The normalized radii become smaller for shorter wave- 
lengths in each object from the mid-lR to near-lR. 

fer to this spatial scale using the spatial wavelength A = A/B, 
which is the reciprocal of the spatial frequency for the configu- 
ratio n, and corresponds to th e spatial resolution of the configura- 
tion dKishimoto et al.l l2009b; see Sect l3.ll ). It is straightforward 
to show that the spatial wavelength A0.5 at which visibihty be- 
comes 0.5 for a Gaussian is related to R1/2 by 

In 2 

Ri/2 A0.5 (for a Gaussian). (2) 

If we define the half-visibility radius Ryos as 

_ ln2 A().g 

^vo.5 = A0.5 - -TT' (3) 

n 4.5 

then this is equal to /?i/2 for a Gaussian, and can generally be re- 
ferred to as a good approximation of the half-light radius model- 
independently. This radius can approximately be deduced by in- 
terpolating observations when observed visibilities cross over 
0.5 at the longest available baseline. We can also say more gener- 
ically that the visibility at a given spatial wavelength A within 
the first lobe gives the approximate fraction of flux originating 
from the region within the radius ~A/4.5. We note again that for 
a Gaussian brightness distribution, all three quantities, namely 
HWHM, half-light radius R1/2 and half-visibility radius 7?vo.5, 
are identical. 

In the top panel of Fig|5] we compare the visibility functions 
of Gaussian and power-law distributions with various Ri/2- The 
curves are shown as a function of spatial frequency in units of 
cycles per R1/2 (or spatial wavelength in units of R\/2), so that 
we can directly compare different R1/2 power-law cases with 
the same Gaussian curve. We can see that power-law curves 
with smaller values of R1/2 become asymptotically more simi- 
lar to the Gaussian/ring curve. The bottom panel shows the ratio 
of Gaussian HWHM, derived from each visibility, to /? 1/2 as a 
function of the same spatial frequency. For extended power-law 



Fig. 7. Normalized half-light radii Ri/2/Rm (red and green for 
13 and 8.5 fim, respectively), plotted against UV luminosity L, 
defined as 6 yLv(5500A) dKishimoto et al.ll2007h . The red and 
green dotted lines are power-law fits to Ri;2/Rm at 13 fim and 8.5 
fim, respectively. The normalized radii become more compact in 
the objects having the central engine with a higher UV luminos- 
ity. Also shown in purple are thin-ring radii at 2.2 ^um normal- 
ized b y Rin for the sample observed with the KI dKishimoto et al.l 
2011), including two overlapping objects. The gray dotted hori- 
zontal line indicates radius - R^a- 



distributions, the dependence of the Gaussian HWHM on the 
baseline length is quite strong, but for steep distributions, the 
dependence on baseline becomes quite small. All curves cross 
the visibility V of 0.5 at spatial wavelength A close to 4.5 R\/2, 
meaning that the half-visibility radius 7?vo.5 closely represents 
the half-light radius R1/2 in all cases. 

The dotted curves in Fig. [3] are the visibility functions for 
the power-law brightness distributions with three dififerent steep- 
nesses, which are quantified by the normalized half-light radius 
^1/2/^111=3, 10, and 30 (corresponding to power-law indices of 
-2.6, -2.0, and -1.5, respectively). These seem to closely describe 
the mid-IR surface brightness distribution over the sample. For 
the more compact structures seen in the higher luminosity ob- 
jects, it is more diflicult to differentiate the power-law from a 
Gaussian/ring for a given measurement uncertainty. However, 
based on relatively resolved cases such as in ES0323-G77, we 
infer that the power-law description with different values of 7? 1/2 
over the sample is adequate at least approximately. We can then 
translate the observed visibility curves on the 2D plane (Fig|3] 
or the middle panels of FigslT]and|2]i into a value of R1/2 at each 
wavelength, and then disscuss the relation between this size in- 
formation and other observables. 

We derived the best-fit normalized half-light radius R\/2/Rm 
for each object by fitting the visibility curve of a power-law 
brightness distribution to the observed visibility data, separately 
for four wavelength bins, namely 8.2-10, 10-11, 11-12, and 12- 
13 /im in the observed frame. To ensure a uniform comparison at 
the same rest wavelength over the sample, we fitted a power law 
as a function of wavelength over the derived Ri/2/Rm (and its er- 
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Table 6. Estimated half-light radii R1/2 in two rest-frame wavelengths in the mid-IR and their uncertainties in dex. 



name 


"ill 




O.J //m 






13 fivn 






(pc) 


RljllR'm R 


1/2 (pc) 


unc. (dex) 


Ri/2/Rm 


Ri/2 (pc) 


unc. (dex) 


NGC4151 


0.033 


15.9 


0.52 


0.18 


36.8 


1.21 


0.14 


NGC3783 


0.061 


10.0 


0.60 


0.09 


46.9 


2.84 


0.20 


ES0323 


0.10 


4.7 


0.49 


0.08 


7.9 


0.83 


0.12 


H0557-385 


0.12 


4.1 


0.51 


0.21 


10.2 


1.26 


0.21 


IRAS09149 


0.47 


2.2 


1.06 


0.09 


2.8 


1.33 


0.13 


IRAS 13349 


0.49 


3.6 


1.76 


0.12 


3.2 


1.58 


0.20 


" SeeTableUJfor^in (= 


Rj^ fit) in mas. 













ror) at the four wavelength bins, and used this fit to derive values 
at 8.5 and 13 jum in the rest frame. The results are summarized 
in Table|6] The uncertainties are formal 1-cr uncertainties for the 
fits, except for NGC4151 and NGC3783 where the targets are 
well-resolved (V < 0.5 at long baselines) and the uncertainty in 
Ri/2 was estimated as the dispersion in those derived from differ- 
ent baseline data. Since the dispersion in /?i/2 about its best- fit 
value was more symmetric in log space, the uncertainties are 
shown in dex. 

We note that the fits were done with the inner and outer ra- 
dius fixed at 1 and 100 Rin- The latter can affect the radius es- 
timation for the extended cases of NGC4151 and NGC3783. 
However, as long as each source is at least approximately un- 
resolved by the single telescope, the half-visibility radius Ryo.s 
is quite well-constrained as can be seen in the middle panels of 
FiglTJj and FiglU? as well as in Fig|3](see Sect l3.2b . Therefore, 
our Ri/2 estimation would also be quite robust for the two tar- 
gets. 



4. Discussion 

4.1. Sizes at various wavelengtiis as a function of UV 
luminosity 

Using the half-light radius R1/2 derived in the previous section, 
we discuss below the overall relations of the size to observing 
wavelengths and flux-related quantities such as spectral shape, 
total flux, and luminosity. 

Fig|6] shows Ri/2/Rm as a function of rest-frame wave- 
lengths. For the near-IR, we have plotted the thin-ring radii 
from the KI data for the two obj ects in our present MIDI sam- 
ple, NGC4151 and IRAS13349 (iKishimoto et al.ll2009ah . The 
thin-ring approximation should be fairly robust in the near-IR, 
where the radial brightness distribution is expected to be much 
steeper and the power-law visibility curve becomes very close 
to that of a Gaussian/ring (Fig|5]l. A more detailed estimation 
of /?i/2 would require knowledge of the shape of the innermost 
dust distribution, which is not yet well constrained. Thus, we 
simply adopt here the thin-ring approximation for the near-IR 
data. We would infer that the other four targets, for which no 
KI data are available, probably have near-IR ring radii of ~l-2 
J?in based on the Kl-ob served sample of type 1 AGNs (see Fig|7J 
IKishimoto et al.ll201 ih . From the mid-IR to the near-IR, the nor- 
malized emission size would then become certainly smaller with 
wavelength, very roughly following some power law for each 
object, but with different indices for different objects. 

As we demonstrated in Figl3] the mid-IR emission size 
seems to have a luminosity dependence. This is now quantita- 
tively shown in Fig|7] where R\/2/Rm is plotted against UV lu- 
minosity L. The data points are color-coded for the observing 
wavelengths, and we include the 2.2 yum thin-ring radii of all the 



Kl-observed sample from 'Kishimoto et al.l (1201 lb . As we can 
see, the ratio of half-light radius to R^n becomes smaller with 
increasing UV luminosity. This implies that, for a given radial 
temperature distribution, the higher luminosity objects have a 
steeper radial density structure, at least in the mid-IR emitting 
radii, as we discuss in greater detail below. We also plot power- 
law fits to the Ri/2/Rm at 8.5 and 13 ;um in Fig|7] 

We can also analyze the luminosity dependence of 
Mnnormalized /?i/2 in a physical scale. Fig|8]shows the half-light 
radius 1/2 in pc plotted against the same UV luminosity. In the 
physical scale, the radius in the mid-IR increases with luminos- 
ity much slower than L'''^. Power-law fits to R1/2 in pc at 8.5 and 
1 3 fim as a function of luminosities are 



/6vLv(5500A)\ 

Ri/2 (8.5//m) = (1 .2 ± 0.2) • 

\ 10™ erg/ sec / 



0.21 ±0.05 



and 



7?i/2(13;um) = (1.3 + 0.3)' 



6vL^(5500A) \ 
1046 erg/sec j 



0.01 ±0.07 



pc, 



(4) 



(5) 



where we see that R1/2 at 13 fum is consistent with being con- 
stant. 

These results are contrary to those of previous studies of 
the m id-IR size as a function of lumino sity by Tristram et al] 
(l2009h and lTristram & SchartmannI (l201ll) . who concluded that 
the mid-IR size is consistent with being proportional to L'^^. 
They used a single Gaussian for the size estimation, and this 
can lead to a large systematic error owing to the baseline depen- 
dency. When the radial brightness distribution has a relatively 
shallow power-law form, a Gaussian size would represent the 
emission size that is resolved by a given interferometer con- 
figuration, leading to the strong dependence on the sampled uv 
points. Thus, in this case, the simple Gaussian size is quite in- 
adequate for representing the emission size for the whole distri- 
bution. Here we account for the effect of the uv sampling more 
properly, based on the inference that the brightness distribution 
is of roughly a power-law form. This enables us to perform a 
more accurate, uniform comparison over the sample, as shown 
in Fig[8] To gain a greater accuracy and model-independency, 
we can always return to the uniform view of the visibility curves 
over the sample shown in Fig|3] 

The half-light radius in units of Ri„ has a quite close relation- 
ship with the mid-IR spectral shape. Figure |9] shows Ri/2/Rin 
plotted against spectral index a in fy oc ■/' in the mid-IR. The 
shape becomes bluer when the size is smaller In terms of the 
radial structure, this is what we would expect, as the steeper ra- 
dial structure looks bluer as the contribution from hotter dust 
becomes more dominant. We elaborate on this point later using 
simple models. 
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Fig. 8. Half-light radii R1/2 in pc (red and green for 13 and 8.5 fim, respectively), plotted against UV luminosity L, which is defined 
as 6 vLv(5500A) (.Kishimoto et al.ll2007l) . The red and green dotted lines are power-law fit to R1/2 at 13//m (oc [ji ^i^^o.oi^^ ^j^^ 8.5jum 
(oc /^o.2i±o.05-)^ respectively. In this physical scale, these mid-IR radii increases with luminosity much slower than L''^, or almost 
constant at 13 fim. Thin-ring radii at 2.2yum also in pc are plotted in purple for the Kl-observed sample. Plotted with gray plus signs 
are the near-lR reverberation radii (Suganuma et al. 2006 and references therein), with the black dotted line showing their L'^^ fit. 
This fit is the definition of the inner radius Ri„ in this paper. 



We can also calculate an average surface brightness in the 
rest frame using R1/2 as {vLy I An) I (nR^^i^, which is equivalent to 

Vobs/v„te • (1 + zfl{n6\i^) where 6112 is the angular size of T^i/a. 
FigfTOl shows this as a function of rest-frame wavelengths. This 
can be directly compared with the Planck function as shown 
in the same figure, where the comparison gives the brightness 
temperature' at a given wavelength. In the near-lR, all the ob- 
jects show the brightness temperature ~1400-1000K, while in 
the mid-lR the brightness temperature ranges from ~1000K to 
~200K. In the near-IR, the brightness temperature is quite close 
to the observed color temperature of ~1400K, meaning that the 
average surface filling factor at the innermost dusty region must 
be close to unity. Figure [TT] shows the surface brightness as a 
function of UV luminosity, and the corresponding brightness 
temperatures at 2.2 and 13 /vm are shown in the second and third 
y-axis, respectively. We can again clearly see the luminosity de- 
pendence. 

The close relationship between Ri/2/Rin and the mid-IR 
spectral index (Fig|9]l and the relationship between Ri/2/Rm and 
UV luminosity (Fig|7]) implies that there should also be a corre- 
lation between the mid-IR spectral index and luminosity. In fact, 
the luminosity de pendence of the radi al dust distribution was al- 
ready inferred by iHonig et"an (1201 Oh . based on the IR spectral 



comparison between Seyfert galaxies and high-luminosity type 
2 QSOs. Here we directly show the distribution change by spa- 
tially resolving it. 

4.2. Multiple-temperature Gaussian fit 

While the power-law brightness description at each wavelength 
gives a relatively model-independent measurement of the emis- 
sion size, it does not physically describe the change in the over- 
all size with wavelength and the flux that we observe. A sim- 
ple, though very approximate way of describing the power-law 
brightness with changing steepness over different wavelenths is 
to consider two or more rings/Gaussians of different size at dif- 
ferent temperatures. At each wavelength, the multiple compo- 
nents with different sizes contribute differently to the visibil- 
ity function, giving a wider range in spatial frequency than the 
range covered by a single ring/Gaussian. The contribution ra- 
tio of the multiple temperature components changes with wave- 
length, causing the overall apparent size to change with wave- 
length. 



' Here the brightness temperature is defined without requiring the 
intensity to be in the Rayleigh- Jeans hmit. 
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Table 7. The results of two-temperature Gaussian fits. 





(pc) 










mid-IR Gaussian 


' TTi inor/tn n lor 


T(K) 


HWHM (i?i„) 


fo 


T(K) 


HWHM (R„,) 


fo 


NGC4151 


0.033 


1389+68 


1.06+0.02 


0.36+0.10 


282+6 


51.1 + 1.2 


0.12+0.01 


12.9 


NGC3783 


0.061 


1220+108 


(fixed to 1) 


0.57+0.28 


242+9 


32.4+0.6 


0.77+0.14 


8.1 


NGC3783" 


0.061 


1220+108 


(fixed to 1) 


0.57+0.28 


242+9 


39.5+1.0 


1.12+0.21 


4.1 0.46+0.02 


ES0323 


0.10 


1539+175 


(fixed to 1) 


0.29+0.14 


320+18 


12.3+0.4 


0.71+0.15 


1.7 


H0557-385 


0.12 


1241+54 


(fixed to 1) 


0.97+0.26 


307+15 


18.2+1.2 


1.06+0.26 


1.4 


IRAS09149 


0.47 


1204+103 


(fixed to 1) 


0.63+0.32 


290+39 


7.4+0.5 


1.67+0.89 


0.65 


IRAS 13349 


0.49 


1314+57 


1.57+0.06 


0.37+0.10 


366+26 


9.6+0.7 


1.9+0.6 


0.98 



" An elliptical Gaussian is fit for the mid-IR component, assuming the major axis at 136°. The quoted mid-IR HWHM is along this major axis. 
* See Table[T]for R,^ (= Rr^ fit) in mas. Reduced values for visibility fits. 
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Fig. 9. The mid-IR (8.5-13 fim) spectral indices in fy are plot- 
ted in red against normalized half-light radii Ri/2/Rm at 13 fim. 
Also plotted in purple are spectral indices between 2.2 and 13 
yum against thin-ring radii at 2.2 fim normalized by Ri„. For the 
four objects without near-IR interferometry data, we indicated 
the near- to mid-IR spectral indices in dashed lines over the thin- 
ring radius range of 1-2 roughly expected based on other 
targets observed with the Keck interferometer (see Fig|7]i. The 
dotted line shows the grid for a simple power-law model with 
radial surface density index a of (-2.0, -1.0, -0.0) and temper- 
ature index /3 of (-0.5, -0.4, -0.3) with the innermost tempeture 
rin=1400K. The grid covers the observed range of Ri/2/Rm, but 
systematically bluer in spectral shape. The dashed lines show 
the same grid but with rin=700K, which apparently matches 
the mid-IR observations over the sample when the density in- 
dex ranges from — 1 to in different objects with a temperature 
index of — 0.35 + 0.05 (for details, see Sectl4!4li. 



We adopted a Gaussian geometry for each component, 
though rings might provide a more physically accurate descrip- 
tion of in particular the innermost dusty region. In terms of the 
visibility function, a ring and Gaussian essentially give the same 
curve at the low spatial frequencies before the first null (FigHJ 
see e.g.|Millour 2008). Both a thin-ring and a Gaussian geometry 
have only one parameter, but a Gaussian permits us to more sim- 
ply quantify the average surface brightness of the component. 



Fig. 10. Average surface brightnesses vly in the rest frame, de- 
fined as {vLy I An) I {jiR\ are plotted against rest wavelengths. 
Overplotted are the Planck functions with temperatures 1400, 
700, and 350 K (dotted, dashed, and dash-dot, respectively). For 
NGC4151 and IRAS13349, which have near-IR Keck interfer- 
ometry data, the average surface brightness at 2.2 //m is cal- 
culated in the same way but using thin-ring radii. For the ob- 
jects without near-IR interferometry, it is calculated assuming 
/?i/2=(1.5+0.5) Rin and indicated by dashed lines. See further 
details in Sectl4TT] 



and its visibility curve (which is also a Gaussian) is the simplest 
(no nulls, no lobes). 

We describe each Gaussian component as a brightness dis- 
tribution 5 V at a radius r from the center given as 



5 v( V, r) = B,(v, r) ■ /o In 2 ■ e"^' ^"'m , 



(6) 



where By is the Planck function at frequency v and temperature 
T, and 7? 1/2 is the physical size of the Gaussian in HWHM, where 
HWHM - R\/2 for Gaussian. The isotropic luminosity Ly of this 
component is 



Ly/An = fo ■ By(v, T) ■ kRI^. 



(7) 



The factor fo gives the average surface brightness with respect to 
the blackbody when all the emission is contained within HWHM 
(see more on emissivity and surface filling factor described in 
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Fig. 11. The same average surface brightness as in Fig[TOl but 
plotted against the UV luminosity of the central engine, defined 
as 6 vLy(5500A) dKishimoto et'aP 120071) . The red and green 
symbols correspond to 13 and 8.5 fim, respectively. The surface 
brightness at 2.2 /im is sh own in purple for t hose w ith the near-IR 
interferometry data from 'K ishimoto et al.l (1201 ih . The luminos- 
ity dependence is quite clear at least for the surface brightness in 
the mid-IR. 



Sect l4.4l ). In addition to these two Gaussians, we also include 
the accretion disk component, which is assumed to remain un- 
resolved and have a spectrum of /y oc v'^^^^ at IR wavelengths 
longward of 0.8 fim ('Kishi moto et'al]|2009al l2008V 

Table |7] shows the results of two-Gaussian-component fits. 
We assumed that HWHM= 17?in for the ~ 1400K near-lR compo- 
nent here, except for the two objects for which we have Kl data. 
We see that the ~300K mid-lR component has a HWHH of sev- 
eral to a few tens of Ri^, which implies an effective temperature 
radial gradient of /? from — 0.8 to ~-0.4 where T oc r^. This 
effective index includes the effect of the radial gradient of sur- 
face density distribution. In section 14.41 we attempt to separate 
the temperature and density gradients using more physical but 
simple models. 



4.3. PA dependence 

The PA coverage of our baselines for each object is very limited. 
Since our targets are all type 1 AGNs, we do not expect to see a 
large PA dependence in visibilities. Nevertheless, we seem to see 
some PA dependence at leas t in one obje ct, NGC3783. This has 
already been pointed out bv lHonig et al.l (12010). The brightness 
distribution along the SE-NW direction seems more extended 
than that in the NE-SW direction. 

We can complement our inner PA dependency study with op- 
tical polarization measurements. The projected direction of the 
system polar axis can be inferred to be either parallel or per- 
pendicular to the PA of the optical continuum polarization, de- 
pending on the line polarization properties ( Smith et al. 2005 
iKishimoto et al.l 120041) . If a linear radio structure is detected, 
this direction also indicates that of the system polar axis. All 
the targets except IRAS09149 have optical polarizatiori mea- 
surements (Brindl e et al1ll990t IWills et al.llT99l lMartellll998l: 
[Schmid et al. 2003^, and a clear radio axis is also available for 
NGC4151 (jMundell et al.l l2003h . The resulting inferred PAs are 
summarized in Table |5] and shown in the bottom panel of the uv 



coverage for each object in Fig [T]and|2] with polar and equatorial 
directions in dashed and solid lines, respectively. 

To evaluate the elongation of the mid-IR emission in 
NGC3783, we attempted to fit an elliptical Gaussian for the mid- 
IR component with the major axis fixed at the polar axis PA 
136°, and the results are listed in Table|7] The emission seems 
to be a factor of ~2 elongated along this polar direction, sug- 
gesting that some of the mid-IR emission originates from the 
polar region. More extensive data and analyses will be presented 
for this object by Honig et al. (201 1 in prep). 

4.4. Power-law brightness distribution: 
temperature/density-gradient model 

4.4.1. Model description 

We now discuss the direct physical constraints on the radial 
structure of the inner dusty region thermally emitting in the IR. 
To this end, we use a simple but generic description of the radial 
surface brightness distribution with as few parameters as pos- 
sible. To physically describe a face-on radial brightness distri- 
bution Sv{v, r) at radius r and observing frequency v, we would 
need at least two separate distributions, namely the temperature 
and surface density distributions. Here we slightly generalize the 
surface brightness description in IKishimoto et al.1 (l2009h) as 



5v(v,r) = /,(r)-B,(v, W)). 



(8) 



The idea behind this description is that, at each radius, the face- 
on surface brightness is dominated by the emission from the 
heated dust g rains that h ave a maximum temperature at that ra- 
dius (e.g. Honig & Kishimo to 2010). These grains are probably 
at around the surface of the torus, and most likely directly illu- 
minated by the central source, if there are such directly illumi- 
nated grains at that radius. Under this maximum-temperature ap- 
proximation, we refer to the dust at Tmax simply as heated dust 
grains/clouds. We note that our observations are insensitive to 
much colder, interior dust grains that may exist deeper along our 
line of sight at that radius.^ We write this maximum temperature 
Tm-ixif) at radius r as 



Tm ■ {rlRinf- 



(9) 



We mainly consider a discrete, clumpy, or inhomogeneous 
cloud distribution. The surface brightness of each heated cloud 
at radius r would approximately scale with r in proportion to 
BviTm-ixir))- The additional factor f^ir) above, which we write 
in the form of a power law as 



feXr) = /o ■ {rlRinf, 



(10) 



describes the radial run of the surface density (i.e. per unit area) 
distribution of these heated discrete clouds (i.e. at T^^^), but 
also includes the emissivity of each heated cloud surface. More 
specifically, the dimensionless factor can be regarded as a sur- 
face filling factor multiplied by the emissivity. In the case of 
optically thick illuminated clouds where their emissivity does 
not depend sensitively on the radial distanc e from the illuminat- 
ing so urce or the observing frequency (e.g. lHonig & Kishimotol 
l20ld, see their Fig. 3), the factor /^(r) is roughly proportional 
to the radial surface density distribution of the heated clouds, 
or heated dust grains. Thus, as long as the emissivity is roughly 



- There should be very little or no directly illuminated dust grains 
situated at greater depth from the surface along our line of sight, because 
the vertical optical thickness (along our line of sight) is quite likely 
smaller than the equatorial optical thickness. 
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Fig. 12. (a) SED in vLy for each object normalized at /I = 2.2 yum, plotted against normalized half-light radius /?i/2//?m, for 2.2, 8.5, 
and 13 yum (purple, green, and red, respectively), (b) The same data but with the grid in dashed lines for temperature/density gradient 
models with radial surface density indices of a=-0.4,-0.8,-1.2,-1.6 for the case of T oc r plotted for the three wavelengths 
2.2, 8.5, and 13 //m (purple, green, and red squares). This simple power-law model does not explain the observed SED and size at 
the same time, (c) The same data but with the grid in dashed lines for the power-law (density oc ^-0-0.2,-0.4,-0.6^ j ^ ^-0.36-) _j. ]^400K 
thin-ring at R\a, which roughly explain the observed SED and size. The two reddened sources IRAS 13349 and H0557-385 could be 
interpreted in terms of the same model with a slightly larger filling factor for the power-law component. 



Table 8. The fit results of temerature/density gradient model + inner ring. 



name 


power-law {T,„ = 700 K) 




inner 1400K ring 


■ 




polar / 
equatorial 


density index a 


T index yS 


,/o 




radius {R\^ 


,/o 




NGC4151 


-0.31 ±0.04 


-0.26 ± 0.01 


0.16 ± 


0.01 


1.25 ±0.03 


0.26 ± 


0.02 


3.5 




NGC3783 


+0.04 ± 0.05 


-0.31 ±0.01 


0.16 ± 


0.02 


(fixed to 1) 


0.29 ± 


0.02 


2.4 




NGC3783" 


+0.04 ± 0.06 


-0.35 ± 0.01 


0.13 ± 


0.02 


(fixed to 1) 


0.28 ± 


0.02 


0.69 


2.43 ±0.18 


ES0323-G77 


-0.77 ± 0.09 


-0.30 ± 0.02 


0.55 ± 


0.09 


(fixed to 1) 


0.42 ± 


0.03 


0.60 




H0557-385 


-0.65 ±0.10 


-0.31 ±0.02 


1.40 ± 


0.26 


(fixed to 1) 


0.47 ± 


0.02 


1.3 




IRAS09149 


-0.46 ± 0.21 


-0.42 ± 0.03 


0.69 ± 


0.19 


(fixed to 1) 


0.27 ± 


0.02 


0.68 




IRAS 13349 


-1.26 ±0.16 


-0.27 ± 0.03 


3.6 ± 


1.0 


1.87 ±0.07 


0.17 ± 


0.01 


0.71 





" Fit with additional ellipticity for power-law component, assuming polar axis at 136°. 



constant over the radius and frequency, we can regard /e(r) as a 
normalized surface density of heated dust grains (i.e. again those 
at Tmax)- In this way, our description is generic enough to ensure 
that we can obtain as direct constraints as possible from the ob- 
servations. 

The model would be called a temperature gradient model, 
if the density distribution is set to be uniform, i.e. a = 0, and 
the brightness described only by the temperature index fi. In this 
case, the index would be an effective temperature index, which 
would include the effect of the density gradient. Here we try to 
separate the temperature and density gradient. At least in the 
case of the direct illumination of a dust grain, we would expect 
that /3--Q.5 for large grains and /3 — 0.36 for standard-size ISM 
grains (e.g. Barvainis 1987). By refeiTing to these values, we 
attempt here to obtain constraints on the surface density distri- 
bution. 



4.4.2. Mid-IR spectral slope and SED-size relation 

In Fig|9l we show a grid for the two indices a and /3 on the plane 
of the mid-IR spectral index and half-light radius. The upper- 
most grid with dotted lines is for Tin = 1400K, which corre- 



sponds to the observed color temperature in the near-IR and is 
thought to be roughly equal to the dust sublimation tempera- 
ture. The range of density index a from -2 to for a given 
temperature index reproduces the observed range of normalized 
half-light radii Ri/2/Rm- However, the expected mid-IR spec- 
tral shape as shown by the dotted grid is systematically bluer 
than observed. In fact, the observed range of the mid-IR spectral 
index, as well as that of Ri/2/Rin, is reproduced well with the 
power-law distribution if Tin is ~ 700K, where the correspond- 
ing grid is shown in dashed lines in Fig|9l We need to accommo- 
date this size-color constraint. 

The total flux spectra, including the mid-IR spectral shape, 
and the size information for all the objects can be summa- 
rized on the plane of the SED, normalized at 2.2 fim where 
all objects show similar brightness temperature (FigfTOt. plot- 
ted against R\/2/Rin as shown in FigfTSb. For objects without 
Keck near-IR interferometric data, we set their R1/2 at 2.2 //m 
to have the expected r ange of values from 1 to 2 Rin (see Fig|7j 
iKishimoto et alJl201 ih . In the mid-IR, the objects show a rela- 
tively uniform behavior. From the high to low luminosity ob- 
jects, as the normalized size R\/2/Rin becomes larger, the mid-IR 
spectral shape becomes gradually redder (which is already seen 
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Fig. 13. The modeled surface brightness distribution of the 
power-law component at 13 fim as a function of radius in pc. 
The value is normalized with the constant value of the Planck 
function at 13 /vm for r=1400K. The color indicates the UV lu- 
minosity of each object as in Figl3l which illustrates that higher 
luminosity objects have steeper and brighter surface brightness. 



in Fig|9]), and the mid-IR total flux level gradually increases with 
respect to the near-lR flux, although two objects, IRAS 13349 
and H0557-385, seem to have slightly (a factor of a few) higher 
mid-lR flux level than the other four objects (see more below). 

Overplotted in dashed lines in FigfTSb on the same data 
set is the SED-size relation at three wavelengths expected from 
the power-law distribution with the temperature index j6=-0.36 
over the range of density indices a=-l.6 to -0.4 (the case with 
slightly steeper /3 is similar, but with a slightly different den- 
sity index range). The total flux mismatch can be clearly seen 
across this plane. It is quite difficult to reproduce the wide range 
of RijilRm and the overall near- to mid-lR SED which is rela- 
tively flat in v/v, while keeping the relatively red mid-lR spectral 
shape, for this single power-law structure. 

4.4.3. Mapping the radial structure witli a briglit near-IR rim 

However, we would be able to roughly understand both the total 
flux and size information if, in addition to the power-law struc- 
ture, there is a central brightness concentration on the scale of 
~Rin emitting nearly at the sublimation temperature of ~1400K. 
This concentration would correspond to a dust sublimation re- 
gion, and could be a bright rim of the innermost dust distribution, 
such as the puffed-up rim that is of ten discussed for the inner re- 
gion of young stellar objects (see iDullemond & Monniejl2010l 
for a review). The exact structure of this central concentration is 
not yet well constrained with our data. FigfT2b shows a grid of 
models with the power-law plus the central concentration repre- 
sented simply by a 1400K ring at r-Rm, which seems to closely 
reproduce the overall SED-size relation with a range of density 
indices. In this case, the two objects, IRAS 13349 and H0557- 
385, simply have a slightly higher flux level in the power-law 
component with respect to the near-lR ring. 



Fig. 14. The modeled surface density distribution of heated dust 
(at maximum temperature at each radius) in the power-law com- 
ponent, illustrating that higher luminosity objects have steeper 
and denser distributions. 



In Figures [T] and |2] we also show the best-fit results with 
this model for each object and Table [8] summarizes the parame- 
ters. The normalization of the ring luminosity is parameterized 
in the same way as in the case of Gaussian given in Eq|7] where 
Ri/2 is replaced by the radius of the ring. Figure [13] shows the 
resulting surface brightness profiles of the models for each ob- 
ject. The power-law component in higher luminosity objects be- 
comes steeper, as expected. Figure [14] shows the corresponding 
normalized surface density distribution of the heated dust grains 
(i.e. with a maximum temperature at each radius). It becomes 
steeper, and also slightly denser, for higher luminosity objects. 
The density indices are from ~0 for lower luminosity to — 1 for 
higher luminosity objects in our sample. 

We note that the central temperature Ti„ of the mid-IR power- 
law distribution has to be as low as ~700 K, to reproduce the 
relatively red mid-IR spectral shape. (We fixed Ti„ to be TOOK 
in the fits implemented above, and we would need visibility data 
that fill the gap between near-IR and mid-IR wavelengths to con- 
strain Tin more tightly.) This would mean that there is a compact 
region of such a low temperature quite close to the radius Rin, 
which could be the shadowed region just behind the bright sub- 
limating rim. 

4.4.4. NGC4151 

In the case of NGC4 1 5 1 , this low-temperature compact region is 
probably what we see in the mid-IR visibility curves, which be- 
come relatively flat at spatial wavelengths smaller than ~ 100 Rin 
as shown in the mid-panel of Fig[T]3. In our interpretation, this 
is caused by both the high spatial frequency part of the low-Tin 
power-law component and the hot central brightness concentra- 
tion. Here the effect of the latter on visibility is indicated in this 
panel (also in the other mid-panels in Fig [T] and [2] for other ob- 
jects), where the dotted curves show the visibilities without the 
bright near-IR rim. 
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These can also be understood in terms of the correlated flux 
observed at long baselines (see the top-panel of FiglTJj), which 
provides an approximate spectrum of the unresolved part of the 
source. Its color temperature is ~400 K (or ~500 K if we exclude 
slightly less reliable data at wavelengths shortward of ~10 yum), 
which is much higher than that of the total flux (~280 K) and 
much lower than the sublimation temperature. 

The deviation of our model fit from the data might indicate 
that the low- Tin core has even more flux while the outer radiat- 
ing source is even more extended. This could suggest a model 
of three compo nents (including the n ear-IR component), which 
was implied bv Burtsc her et al.l (12009). who modeled their two- 
baseline mid-IR data with a Gaussian and an unresolved source, 
and suggested another near-IR component. We note, however, 
that, puzzlingly, the correlated flux in their two-baseline data did 
not show the relatively blue color temperature, thus gave almost 
no indication of the size becoming smaller toward shorter wave- 
lengths within the mid-IR. 

4.4.5. Warm mid-IR and hot near-IR parts 

The composite brightness structure discussed above, consisting 
of emission components of a warm mid-IR dominating part and 
a hot near-IR dominating part, also seems to be discernible in 
the SED, particularly in those of relatively red objects such as 
NGC3783 and NGC4151 (FigHJ. The two-component explana- 
tion of the near to mid-IR SEDs has been advoc ated by different 
authors from SED analyses (e.g. lMor et al ]|M)9). Historically, 
the bump around 3 yum in the SEDs of QSOs, which could corre- 
spond to the hot near-IR part, has be en known for decades (e.g. 
see sec.V in lNeugebauer et al ] ll979h . 

A direct size constraint on the hot central brightness con- 
centration can be derived from near-IR interferometry, and its 
structure might be a little distinct from the warm power-law dis- 
tribution. For example, the near-IR visibility of NGC4151 sug- 
gests that it has a quite compact structure with R\i2 quite close to 
Rin, while most of its mid-IR emission is quite extended. In that 
case, another parameter, in addition to the heating luminosity, 
may be responsible for forming the structure of the inner region. 
We plan to explore these regions with further measurements in 
both the near-IR and mid-IR. 



5. Conclusions 

We have presented mid-IR long-baseline interferometry of six 
type 1 AGNs. Within the small sample, which nevertheless spans 
over ~2.5 orders of magnitudes in the UV luminosity L of the 
central engine, we have shown that the radial structure of the dust 
distribution changes systematically with luminosity. We have ar- 
gued that the mid-IR radial brightness distribution can approx- 
imately be described by a power-law with the inner boundary 
set by the dust sublimation radius Ri^. In this case, the power- 
law steepness can be quantified with a half-light radius /? 1/2. The 
mid-IR half-light radius ^1/2 in units of becomes smaller (i.e. 
the power-law distribution becomes steeper) as the luminosity 
increases, with /? 1 /2 at 1 3/^m ranging from a few tens of Ri^ down 
to a few Ria over the range of L from ~ lO^-''^ to ~ lO^*" erg/sec. 
This means that, in contrast to the results of previous studies, 
the physical mid-IR radii in pc is not proportional to L'''^, but 
increases with L much more slowly. Our current estimate is that 
the radii is oc /^+o.2i±o.o5 ^j. nearly constant at 13 im\. 

This mid-IR size normalized by is correlated with the 
mid-IR spectral shape, where the normalized size is smaller for 



bluer spectral shapes. The mid- to near-IR SED as well as the 
measured size in mid- and near-IR wavelengths can roughly be 
explained in a picture where there is a central brightness con- 
centration at the scale of ~Ria emitting nearly at the dust subli- 
mation temperature of ~1400K, over an underlying power-law- 
like structure extending outwards. The latter structure is inferred 
to have an innermost temperature as low as ~700K, and a ra- 
dial surface density distribution of the heated dust that becomes 
steeper from lower to higher luminosity objects, ranging from 

to ~r"' . The radial temperature run is consistent with T oc 
of P 0.35 + 0.05. The exact structure of the innermost con- 
centration will further be constrained with the planned near-IR 
interferometry. 
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Appendix A: Reduction and calibration of IVIIDI data 
for faint targets 

We used a part of the software EWS ( lJaffell2004t version 1.7.1) 
to reduce the MIDI data described in this paper. The mid-IR cor- 
related flux of all the targets here is < 0.5 Jy, which is much 
fainter than normally handled with the software. Therefore, we 
implemented several additional procedures to reduce and cali- 
brate the data as described below, using our own IDL codes. In 
short, we simply averaged or smoothed over a relatively large 
number of frames, designated as w here (typically w ~20-40) to 
determine both the group-delay and phase-off'set tracks, and ap- 
plied the same averaging to the calibrator frames to compensate 
for the side efi'ects of the averaging. 

We also observed several calibration stars, listed in 
Table I A. fl which are fainter than ~0.8 Jy at 12 /im, have known 
spectral type, and are unresolved with VLTI/MIDI. The deduced 
correlated flux is directly compared with th e total flux, which 
is obtained by scaling the template spectra of Cohen et al.|[T999l 
(not from less accurate total flux frames obtained with MIDI). 
This facilitates a reliable evaluation of the accuracy of our re- 
duction and calibration procedures. The scaling factors are de- 
rived from both IRAS and AKARI measurements. On the basis 
of the close match between the two, we estimate that the total 
flux accuracy is ~A%. 

A.1. Large gsmooth parameter 

When the group delay is determined from the delay function (the 
Fourier transform of each fringe spectrum from the frequency 



domain to the delay domain), several frames are averaged or 
frames are smoothed over time direction as specified by gsmooth 
parameter in EWS to suppress instrumental delay peaks. The de- 
fault is four, but we used a larger number (typically 10-20) to 
increase S/N in delay peak determinations, as normally done for 
faint targets. 

A.2. Significant averaging over time for pliase-offset 
determination 

Fig lA.lf a) compares the calibrated correlated flux of the unre- 
solved star HD107485 with the known total flux spectrum (i.e. 
the scaled spectrum described above), and Fig lA.U b) shows the 
deduced visibility (using the known total flux). The raw count 
spectrum of the correlated flux is shown in Fig lA.lf c). The star 
has a peak count rate of ~800 cts/sec in the PRJSM-mode spec- 
trum in this particular observation. In these figures, the dotted 
line shows the spectrum from a default reduction (though with 
the same large gsmooth parameter), which determines phase off- 
sets with a small averaging width. In this case, the spectrum suf- 
fers a red bias, i.e. extra counts at the longer wavelength side. 

This seems to be at least partly due to the determination 
of phase offsets at too low S/N. This can be visualized on the 
complex plane of the broad-band correlated flux counts inte- 
grated over the N-band of ~7 to 13 fim. The phase offsets are 
determined from this integrated broad-band flux. When a large 
enough number of frames are averaged, the distribution of the 
broad-band correlated flux measurements on the complex plane 
shows a (partial) circle as shown in Fig lA.U d). where its radius 
represents the amplitude of the correlated flux with some phase 
offset slowly changing over time. 

A noisy phase-offset determination (and subsequent rota- 
tion using the erroneous phase) would lead to a positive bias 
in the correlated flux (e.g. it would detect non-zero correlated 
flux even if the real correlated flux is zero). The red bias is prob- 
ably caused, since the MIDI count spectra usually have a peak 
at shorter wavelengths, by the phase-offset determination hav- 
ing a larger weight at shorter wavelengths (or more accurately 
count peak), thus leading to the longer wavelength side having 
a stronger positive bias. We note also that, in general, correlated 
flux spectra would also tend to be redder if larger coherence loss 
at shorter wavelengths is not correctly calibrated. 

Therefore, we chose a large w for the phase offset track deter- 
mination, and a boxcar smoothing was applied to the group delay 
track for the same number of frames. By applying the same aver- 
aging to the delay and phase tracks of calibrator frames, we esti- 
mated the low effective system visibility arising from this aver- 
aging process, and calibrated it out from the target frames. As we 
show below, we chose w to maximize the S/N of the phase-offset 
determination and the result is insensitive to the exact value of 
the chosen w. 

A.3. Optimum averaging widtli and evaluation of atmosplieric 
conditions 

Here we quantify the effect of this averaging process as a func- 
tion of w, and describe how we can choose an optimal value 
of w. 

A set of broad-band correlated flux measurements for a given 
w as shown in the complex plane in Fig lA.U d) gives a distribu- 
tion of the correlated flux amplitude measurements. The corre- 
sponding histogram is shown in Fig lA.lf e). We can evaluate the 
distribution as a function of w by calculating a fractional disper- 
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sion cr/m (standard deviation cr divided by mean m) for each w. 
This statistical property would be valid as long as the effective 
number of correlated flux measurements is not low. Here, the 
effective number of correlated flux measurements is the num- 
ber of non-rejected frames (see section IA.5I below) divided by 
w. The reciprocal m/cr corresponds to the S/N of the correlated 
flux measurements per smoothed frame. 

Fig lA.lf g) shows the change of cr/m as a function of w (or 
corresponding averaging time interval, 'accumulation time'; see 
upper axis) for the target and its visibility calibrator The curve is 
cut when the effective number of measurements becomes lower 
than 10. In the small w regime for the faint target, the dispersion 
naturally decreases with larger w as photon noise decreases. The 
photon noise here comprises a background Poisson noise and 
fluctuation due to background subtraction residuals. This can be 
estimated from the noise spectra (a few hundred frames with 
large OPD offsets, taken at the start of each fringe track with 
MIDI) and the fractional fluctuation of the delay function peaks. 
This is shown in gray squares in Fig lA.lf g). and closely matches 
the cr/m curve at small w range. The fractional dispersion cr/m 
starts to increase at some width where the effect of a large time- 
averaging width to smear and effectively reduce the correlated 
flux becomes more dominant over the photon noise decrease. 

For the bright calibrator, cr/m stays relatively constant at 
small w where flux fluctuation slowly decreases (indicated by 
gray plus signs), while the smearing effect of averaging slowly 
increases. The dash-dotted curve is the relative correlated flux 
counts of the bright calibrator (normalized at the smallest w), 
showing a relative effective system visibility at each w. For larger 
w, cr/m increases quickly (system visibility declines quickly) 
when the blurring effect of the large-w averaging starts to domi- 
nate. This gives an indication of the coherence time in the mid- 
IR at the time of the observation (~1.5 seconds in the case of 
Fig lA.lf g)). In the large w regime, the fractional dispersion cr/m 
is thus determined predominantly by the atmosphere at the time 
of the observation, rather than the observed source. 

The coherence time and thus the blurring effect is approxi- 
mately the same between target and calibrator frames when the 
two observations take place close in time and on the sky. In this 
case, we can actually calculate the expected cr/m for the target 
frames as a function of w from the calibrator frames. We sub- 
tract the photon noise (gray plus signs) in quadrature from the 
cr/m of the calibrator (solid plus signs with dashed lines) to ob- 
tain the dispersion caused by averaging, and add it in quadrature 
to the photon noise in target frames. This expected cr/m curve 
for the target is shown as triangle symbols with a dotted line 
in Fig lA.lf g). and closely matches the observed cr/m, meaning 
that atmospheric conditions were very similar. In this case, we 
simply apply the same averaging width for the calibrator frame 
to obtain an appropriate estimate of the effective system visibil- 
ity spectra, or transfer function over wavelengths, for the target 
frames processed with this particular averaging width. 

In Fig lA.lf g). the green circles show the relative calibrated 
broad-band correlated flux (normalized at w giving minimum 
cr/m, or maximum S/N) as a function of w, and the green crosses 
show the relative power-law index of the calibrated correlated 
flux spectrum. The calibrated correlated flux is systematically 
higher, and often redder, at small w because of the positive bias 
from low S/N determination of phase offsets. We have chosen w 
to yield maximum S/N for the target. If w is chosen to be larger 
than this value, the result is quite insensitive to the exact value 
of w, although the effective S/N of correlated flux measurements 
is lower at larger w owing to the smearing effect of averaging. 



A.4. Adjustment for the atmospheric difference 

Fig |A.2l shows a similar case to that shown in Fig lA.ll where 
the count rate is slightly lower and the coherence time is slightly 
shorter The maximum S/N is thus slightly lower However, the 
atmospheric conditions between target and calibrator observa- 
tions seem to be well matched. 

In Figs lA.3l and lA.4l we show the cases where there seems 
to be a difference in atmospheric conditions between target and 
calibrator fringe tracks. In these cases, we try to evaluate the 
difference by comparing the expected cr/m curve for the tar- 
get (triangles with dotted line in Fig |A.3f i) and IA.4f i)). as de- 
rived from the calibrator frames, with the observed curve for the 
target (squares with solid lines). We find that we can fit these 
two curves if we simply scale the time length of the calibrator 
frames (a factor of ~0.7 in Fig |A.3f g). and ~2.5 in Fig |A.4f g)). 
The scaling factor probably corresponds to the difference in the 
coherence time between target and calibrator frames. Therefore, 
we infer that we can obtain appropriate calibrations of the tar- 
get frames if we use this time scaling factor. This correction also 
stabilizes the calibrated correlated flux over different w. 

Fig lA.51 shows the calibrated visibility of the 12 low and 
intermediate S/N observations of unresolved stars as a func- 
tion of the maximum S/N per smoothed frame. We calculated 
the visibility in three different wavelength bins, avoiding the 
range affected by atmospheric ozone features. We confirm that 
the calibration procedure above correctly calibrates the data, 
although we see a slight tendency that, toward shorter wave- 
lengths, the correlated flux is underestimated as the maximum 
S/N decreases. Based on the mean deviation and dispersion of 
these calibrated visibilities, we assign a systematic error in the 
correlated flux as 5% of the total flux, when the maximum S/N 
is at least larger than ~2.5. Since the cases at lower S/N remain 
untested, we did not use the target datasets under this level in 
this paper 

A.5. Refinement of group delay determination 

As described in section lA.ll we use a large gsmooth param- 
eter for faint targets to smooth over a large number of delay 
function frames and thus increase the S/N for delay determi- 
nation and also strongly suppress the instrumental delay peaks. 
However, this also suppresses the atmospheric delay peaks at the 
time intervals where the delay is relatively quickly changing. As 
a result, the delay function image becomes quite 'blobby', and 
strong delay peaks are left only where the atmospheric delay did 
not change rapidly. 

To partially compensate for this effect, we implemented a 
simple iteration to determine the delay track. Namely, we first 
interpolated (over time) between the strong delay peaks to de- 
rive an approximate delay track. Then we used this approximate 
atmospheric delay track (plus instrumental delay) to rotate the 
original fringe spectra in the complex plane, and re-derived the 
delay function smoothed with the same gsmooth parameter This 
makes the resulting delay peaks line up straighter over time, and 
thus recovers a larger number of frames with strong delay peaks. 
We cycled over this process a few times. For calibrator frames, 
this resulted in almost all the frames having strong peaks (i.e. 
even with the large gsmooth parameter). For target frames, sig- 
nificant fraction of frames is recovered, resulting in a good rep- 
resentation of the delay track even in a relatively low S/N case 
as shown in Fig |A.3l 

In low S/N cases, we flagged frames that show delay peaks 
at positions that deviate significantly from the overall track ob- 
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Table A.l. Observation log for sub-Jy unresolved stars. 
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Fig. A.5. Comparison of deduced visibilities for sub-Jy calibra- 
tors. 
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Appendix C: Radial structure plot in linear scale 

In Fig|3]we showed observed mid-IR visibilities of all the targets 
in a single plot to uniformly compare the radial structure of the 
objects in normalized units with spatial frequency in log scale. 
Here in Fig lC.ll we show exactly the same figure but with spatial 
frequency in linear scale, again to facilitate comparisons with the 
data in the literature. 



tained above. This was done by taking the histogram of the 
difference between the determined group delay track with that 
smoothed over many frames (typically ~20-40 frames, set as 
2 X gmoooth). The distribution can roughly be descibed by a 
Gaussian, and frames with large deviations (> 10 cr) were ex- 
cluded. In this rejection process, we did not select frames based 
on the strength of the delay peaks in order to avoid possible bi- 
asing. 



A.6. A note on total flux frames 

For the targets observed with UTs that have total flux of less 
than a few Jy, it is better to implement an additional background 
subtraction using the sky regions very close to the target posi- 
tion. The residual of the primary background subtraction using 
chopped frames can still dominate the target flux for these rel- 
atively faint targets. An optional way to further reduce the sys- 
tematic error from the background subtraction is to iteratively 
reject frames with the sky region counts having large deviations 
from the average. This was implemented in a small number of 
cases when the time fluctuation over the photmetry frames was 
significantly reduced after the frame rejection. With this addi- 
tional sky subtraction, we averaged many sets of data, even from 
different nights, to obtain reliable total flux spectra, as confirmed 
with our VISIR photometry and spectrum. 



Appendix B: Mid-IR data plotted in conventional 
formats 

In Fig[T]and|2]we show our flux and visibility data with normal- 
izations using the inner radius Rin. We show the same data in the 
mid-IR (after reddening corrections; see Sect l2.5b using conven- 
tional units in Fig IB. II and IB. 21 to facilitate direct comparisons 
with the data in the literature. 
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Fig. A.l. The results for the unresolved star HD107485 taken on 2009-05-1 1. (a) Our results for the calibrated correlated flux and 
its error are shown in solid and dot-dash line, respectively, while the correlated flux from a default reduction with EWS 1.7.1 (see 
text) is shown in dotted line. The dashed line is the total flux for this star which is obtained by scaling a template spectrum based 
on the known spectral type and IRAS flux, (b) The same as a but deduced visibilities are shown instead of correlated flux, (c) The 
same as a but raw counts spectra are shown, (d) Broad-band correlated flux measurements integrated over 7-13 jum (smoothed over 
a particular time width w) are shown on the complex plane, (e) Histogram of the amplitude of the broad-band correlated flux from 
the panel d is shown. (/) The phase of the correlated flux as a function of frame number, i.e. time. The same track is repeatedly 
shown with 360° offsets to clearly indicate the track as a function of time. The red points show excluded frames, (g) The ratio of 
the standard deviation cr to mean m of the correlated flux amplitude measurements is shown as a function of smoothing width w 
in frames (lower axis) and in seconds (upper axis). The solid line with squares is for the target (HD107485) while the dashed line 
with plus signs is for the calibrator (see Table lA.ll) . The gray squares and gray plus signs (without any connecting lines) indicate the 
photon noise (see Sect lA.3b estimated for the target and calibrator, respectively. The dotted line with triangles shows the expected 
fractional dispersion cr/m curve for the target, which is estimated from the cr/m curve for the calibrator (i.e. the dahed line with plus 
signs) and the photon noise estimations (i.e. gray squares and plus signs with no connecting lines). The dash-dotted line indicates 
the relative correlated flux counts of the calibrator as a function of w, showing the relative effective system visibility at each w. 
The green circles show the relative calibrated broad-band correlated flux as a function of w. The green crosses show the relative 
power-law index of the calibrated correlated flux spectrum also as a function of w. ih) The delay function on the plane of delay 
versus frame number is drawn in gray scale image, showing the group delay track before the iteration described in Sect lA.51 TB4 
green and red points indicate non-rejected and rejected fringe peaks, respectively, after the iteration. 



HD1 07485 



Fringe SNR Uoximlzotion 
EWS 1.7.1 
Totol flux 




9 10 11 12 

Wavelength (iJ.ru) 



13 



Kishimoto et al.: Radial structure of AGN tori 

HD1 07485 



2.0 



V, 1-0 



Fringe SNR Uoximlzotion..:! 
EWS 1.7.1 




7 8 9 10 11 12 13 
Wavelength (fxm) 




9 10 11 12 13 
Wavelength (/Ltm) 



(a) Correlated and total flux 



(b) Visibility 



(c) Raw counts spectrum 




-6x10^-4x10-2x10'' 2x10" 4x10' 6x10' 
8-13 jtim corr. flux reol (counts/sec) 



E 

S 10 



1 



2x10' 4x10' 6x10' 

-13 correlated flux (counts/sec) 



600 
400 



2000 4000 6000 8000 
Frame nunnber 



(d) Complex broad-band correlated flux 



(e) Histogram of broad-band amplitude 



(f) Phase offset track 



Accumulation time (seconds) 
0.1 1.0 10.0 



I 0.5 



E 
cn 

K 0.0 



md2009051 0_01 h 14m_hdiq7485_fc 
md2009051 0_00h26m_hd94264 




10 100 
Time overoging width (fromes) 



2.0 



(g) sigma/mean ratio as a function of width 



Time (seconds) 
60 80 




I . 



4000 
Frome number 



(h) Delay function and group delay track 



I - 



Fig. A.2. The same as FiglAH but for the data taken on 2009-05-10. 
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Fig. A.3. The same as Fig lA.ll but for the data taken on 2009-05-12. The panel (g) shows the cr/m curve after the correction for t^^ 
time scale difference between the target and calibrator, while the panel (/) shows the curves before the correction. 
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^|g. A.4. The same as FigjA^lbut for another unresolved star HDl 10392 taken on 2009-05-09. 
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(a) NGC4151 (b) NGC3783 (c) ES0323-G77 

Fig. B.l. All the mid-IR data shown in Figs [1] and |2] and summarized in Table |5] are plotted here in more conventional ways for 
(a) NGC4151, (b) NGC3783, and (c) ES0323-G77. Top: the total flux and correlated flux spectra taken with MIDI are shown in 
units of Jy. The filled circles indicate the spectra with a 5-pixel boxcar smoothing (3 pixel for NGC4151), while the spectra before 
the smoothing are shown with thin lines. The total flux is shown in black, while the correlated flux is shown in different colors 
for different baselines as indicated by the legend at top-left where the projected baseline lengths and PAs are written. The error 
spectra estimated for the smoothed spectra are shown at the bottom with dash-dot lines using the same corresponding colors. The 
two fluxes from the broad-band VISIR imaging are shown with error bars in thick gray lines. The thin gray spectra mostly at the 
highest flux levels are from the Spitzer IRS observations. The slightly darker gray spectrum shown for ES0323-G77 is from the 
VISIR spectroscopy. For NGC3783, the total flux taken in 2005 scaled to match the flux in 2009 (see Sect j2.1l ) is shown in gray 
color. Bottom: the resulting visibility spectra are shown using the same color scheme, with the error spectra excluding the total flux 
error contribution. The fractional total flux errors are drawn separately in black dash-dot lines. 
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Fig. B.2. The same as Fig lB.ll but the data shown in solid Unes are with the total flux smoothed by 3 pixels while those in filled 
circles are with total and correlated fluxes smoothed by 7 pixels. The figures are for (a) H0557-385, (b) IRAS09 149-6206, and (c) 
IRAS 13349+2438. 
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Fig. C.l. The same figure as Fig[3] but with flie spatial frequency on a hnear scale. 
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